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D ■ ABSTRACT 

Luminous blue variables (LBVs) are a class of highly unstable stars that have been proposed to play a critical role in massive 
stellar evolution as well as being the progenitors of some of the most luminous supernovae known. However the physical processes 
underlying their characteristic instabilities are currently unknown. In order to provide observational constraints on this behaviour 
we have initiated a pilot study of the population of (candidate) LBVs in the Local Group galaxy M33. To accomplish this we have 
obtained new spectra of 18 examples within M33. These provide a baseline of >4 yr with respect to previous observations, which 
is well suited to identifying LBV outbursts. We also employed existing multi-epoch optical and mid-IR surveys of M33 to further 
constrain the variability of the sample and search for the presence of dusty ejecta. Combining the datasets reveals that spectroscopic 
and photometric variability appears common, although in the majority of cases further observations will be needed to distinguish 
between an origin for this behavour in short lived stochastic wind structure and low level photospheric pulsations or coherent long term 
LBV excursions. Of the known LBVs we report a hitherto unidentified excursion of M33 Var C between 2001-5, while the transition of 
the WNLh star B517 to a cooler B supergiant phase between 1993-2010 implies an LBV classification. Proof-of-concept quantitative 
model atmosphere analysis is provided for Romano's star; the resultant stellar parameters being consistent with the finding that the 
LBV excursions of this star are accompanied by changes in bolometric luminosity. The combination of temperature and luminosity 
of two stars, the BHG [HS80] 110A and the cool hypergiant B324, appear to be in violation of the empirical Humphreys-Davidson 
limit. Mid-IR observations demonstrate that a number of candidates appear associated with hot circumstellar dust, although no objects 
as extreme as rj Car are identified. The combined dataset suggests that the criteria employed to identify candidate LBVs results in a 
heterogeneous sample, also containing stars demonstrating the B[e] phenomenon. Of these, a subset of optically faint, low luminosity 
stars associated with hot dust are of particular interest since they appear similar to the likely progenitor of SN 2008S and the 2008 
NGC300 transient (albeit suffering less intrinsic extinction). The results of such a multiwavelength observational approach, employing 
multiplexing spectrographs and supplemented with quantitative model atmosphere analysis, appears to show considerable promise in 
both identifying and characterising the physical properties of LBVs as well as other short lived phases of massive stellar evolution. 
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CN| 1- Introduction agents - specifically binary interaction and impulsive mass loss 

events - may also contribute to this process. 
The advent of recent large transient surveys and deep, multi- Regarding the former, extensive radial velocity surveys have 
. ^ , wavelength observations of the stellar content of external galax- reve aled both a high binary fraction amongst massive stars (e.g. 
^ ■ ies raises the possibility of opening up the final stages of massive Sana & Evans |2QTl1 Ritchie et al. in prep.) and evidence for 
^ ; stellar evolution to unbiased, quantitative study. Since such ob- distinct binary evolutionary channels (e.g. Clark et al. l20Tlbl) 
CO . jects both drive the evolution of their host galaxies while serving cormrrmn g the predictions of previous theoretical studies (e.g. 
as the progenitors of a number of energetic phenomena, includ- p et rovic et al. 121)051 The latter suggestion implicates the mass 
ing core collapse supernovae (SNe), gamma ray bursts and X- loss e pi so des associated with the transitional luminous blue vari- 
and y-ray binaries, an understanding of their brief and violent ables ( LBVs) in mediating the transition to WRs (Smith & 
lifecycle is of wide interest. While it had been supposed that Owocki|2006]l. Observations of Galactic LBVs demonstrate that 
mass loss via (radiatively driven) stellar winds drove the evolu- a subset m assoc iated with massi ve dusty circumstellar nebulae 
tion from the main sequence through the H-depleted Wolf Rayet (_ .01-15 M ; Clark et al. 120031) which have historically been 
(WR), recent observational and theoretical developments sug- tho ught to arise in the so-called 'giant' eruptions; events char- 
ges! that this in fact might be an oversimplification and that other acter ised by increases in bolometric luminosity (L bol ) and which 

help distinguish LBVs from other luminous evolved stars, (e.g. 



* Based on observations made with the William Herschel Telescope Humphreys & Davidson [1994|l. However, recent observations 

operated on the island of La Palma by the Isaac Newton Group in the demonstrate that such outbursts appear to be strikingly diverse 

Spanish Observatorio del Roque de los Muchachos of the Instituto de and are no t alwa ys accompanied by significant mass loss (e.g. 

Astrofisica de Canarias. Clark et al. 2009 ). 
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Given the possible role played by LBVs in the evolution of 
massive stars, there is considerable interest in observationally 
constraining the nature of the outbursts and hence the underlying 
physics driving these instabilities. Unfortunately, the LBV phase 
appears to be brief, resulting in only a handful of widely recog- 
nised Galactic examples (Clark et al. 2005a). However, recent 
IR galactic plane surveys have greatly expanded on the num- 
ber of known candidates (Clark et al. 2003 Gvaramadze et al. 
[20101 Wachter et al. l2010l Stringfellow et al. l20TTb . Such objects 
are invaluable since their spatially resolved circumstellar ejecta 
comprise a record of their mass loss history, although their dis- 
tribution throughout the plane means that long term monitoring 
to determine their duty cycles is an observationally expensive 
proposition. 

An alternative approach is to observe candidates in external 
galaxies. While such a strategy complicates the spatial separa- 
tion of stellar and nebular emission, it is observationally efficient 
given current 4-8 m class telescopes and associated multiplexing 
spectrographs. Moreover, such an approach, informed by ongo- 
ing SNe monitoring campaigns has the potential to efficiently 
identify new examples of giant LBV eruptions for detailed fol- 
low ups; the so-called 'SNe imposters'. In order to complement 
our ongoing program of observations of Galactic LBVs, we have 
recently undertaken a pilot study of (candidate) LBVs in M33 to 
assess the efficacy of this strategy. 

To accomplish this we obtained new 4000-7000 A spec- 
troscopic observations of 17 (candidate) LBVs (Table 1) se- 
lected from the seminal paper by Massey et al. (2007, hence- 
forth Ma07); the resultant >4 yr baseline between observations 
well matching the characterstic timescale of LBV excursions. 
The S/N, spectral resolution and wavelength coverage of the re- 
sultant spectra - encompassing Ha, a prime mass loss diagnos- 
tic, as well as temperature and metallicity diagnostics between 
4000-5000 A - in principle permits non-local thermal equilib- 
rium (non-LTE) model atmosphere analysis of the stars in ques- 
tion. We were also able to supplement this dataset with three ad- 
ditional unpublished spectra obtained in 2003 at the W. M. Keck 
observatory and which are further described in Sect. 2. Of these 
spectra two were of targets already included in our WHT target 
list and one - M33 Var B - was new. 

Moreover, we may also make use of the wealth of existing 
data on the stellar population of M33 in order to investigate both 
the photometric history of, and the presence (or otherwise) and 
properties of dusty circumstellar ejecta associated with the (can- 
didate) LBVs. For this facet of the program we utilised the full 
census of such stars published by Ma07 (their Table 18, not re- 
produced for brevity, which contains all 18 objects for which we 
obtained spectroscopy), supplemented with the remaining WN9- 
1 lh and B super-/hyper-giants within M33 as well as additional 
cLBVs from a variety of sources; the additional stars and asso- 
ciated references are listed in Table J]. 

Consequently, the paper is structured as follows. Data acqui- 
sition and reduction is described in Sect. 2, with the resultant 
spectra discussed in Sect. 3, with the full dataset present online 
in appendix A. Photometric data are presented in Sect. 4 and 
we discuss the nature of the candidate LBVs in Sect. 5, before 
summarising and placing our results in a wider context in Sect. 
6. 
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1 In the absence of a confirmatory spectrum we currently exclude the 
candidate reported by Burggraf & Weis ( 12008 1 . 
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Fig. 1. Multiepoch optical spectra of selected stars Top 

panel - illustrative Iron stars and bottom panel - the P 

Cygni LBV J013416.07+3-3642.1 and the B hypergiant 
J013416.44+303120.8 



2. Observations and data reduction 

The spectroscopic observations of 17 of our targets (Table 
1) were performed at the 4.2 m William Herschel Telescope 
(WHT) in the "Observatorio del Roque de los Muchachos" in 
La Palma, Spain. We used the AutoFib2 plus Wide Field Fibre 
Optical Spectrograph (AF2-WYFFOS) combination for multi- 
object spectroscopy. WYFFOS can be fed with up to 150 1.6"- 
diameter fibres, positioned on the user provided coordinates by 
the AF2 robot over a wide (~ldeg diameter) field of view. 

The observing run started on 2010 September 29 and lasted 
four nights. We designed two AF2 fibre configurations to cover 
our targets. In both cases, we first optimized the number of al- 
located fibers and their priority (which was assigned based on 
obtaining a full sample of each sub-category of candidate LBV 
defined by Ma07, prior indication of variability and optical mag- 
nitude) and maximized the number of allocated sky fibers (see 
below) afterwards. We used the R1200B and R1200R gratings 
to cover the optical-blue spectra (~ 4000-5000 A) and Ha re- 
gion with a resolution of roughly 2 A. The total exposure time 
for the first configuration amounted to 9 hours and 6.5 hours 
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Table 2. Summary of additional candidate LBVs in M33 
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Fig. 2. Multiepoch optical spectra of selected LBVs. Top panel - 
the P Cygni star J013339.52+304540.5; note that the 2010 spec- 
trum is reproduced in Fig. A.3, where weak Sim, On and Cm 
absorption features are visible. Bottom panel - the bona fide LBV 
M33 Var C, with the spectrum of the cool F0-5Ia + hypergiant 
B324 overlayed (Monteverde et al. IT996l and Sect. 5.3.2). 



with the R1200B and R1200R gratings respectively, broken into 
1800s blocks for cosmic -ray subtraction. Due to weather con- 
ditions, the second configuration was only observed during 2.5 
hours with R1200B; consequently, the spectra of stars observed 
in this configuration have much poorer signal-to-noise ratios. 

The data were reduced following the recipe implemented in 
the pipeline provided by the AF2 webpagqj at that time. We 
modified the original procedure to include an initial normaliza- 
tion of the fiat-field before its correction, and an algorithm for 
cosmic ray rejection (van Dokkum 1200 11 1. The code then fol- 
lows the standard steps for multi-fibre spectrum reduction, us- 
ing IRAF0 tasks: bias, fiat-field correction, throughput calibra- 
tion between fibres, wavelength calibration and sky subtraction. 



2 http://www.ing.iac.es/astronomy/instruments/af2/reduction.html 

3 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



Name 


Reference 


M33 Var 2(=VHK2) 


van den Bergh et al. 111975b 


N93351 


Valeevetal. (2009) 


N45901 


Valeev et al. f20T0t 


N125093 


Valeevetal. (2010) 


UIT005 


Urbaneja et al. (120111) 


IFM-B 0301 


Monteverde et al. (11996b 


IFM-B 0515 


Monteverde et al. (119961) 


IFM-B 1040 


Monteverde et al. (11996b 


IFM-B 1054 


Monteverde et al. (119961) 


IFM-B 1330 


Monteverde et al. (11996b 


IFM-B 1345 


Monteverde et al. d 1996b 


OB 10-10 (=UIT 181) 


Massey et al. ( 1995) 



LGGS J013354.85+303222.8 
LGGS J013418.37+303837.0 
LGGS JO 13432.24+304702.7 
LGGS J013432.50+304703.5 



Neugent & Massey 
Neugent & Massey ( 201 1 1 
Neugent & Massey ( 201 1 1 
Neugent & Massey ( 1201 lb 



Top panel: supplemental (candidate) LBVs in addition to those in the 
census of Ma07 (their Table 18) noting that Var83 & B416 are not 
labeled as such in that work, but are LGGS J013410.93+303437.6 and 
LGGS J013406.63+304147.8 respectively. Of these, N45901, N93351 
and N 125093 may be classified as iron stars (Sect. 3) while M33 Var 
2 has historically been classified as an LBV. Middle panel: BA super- 
and hypergiants and bottom panel: Additional Ofpe/WNL and WN9-1 1 
stars from Neugent & Massey J201 II note that these authors now 
classify J013416.07+303642.1 as Ofpe/WNL rather than P Cygni-type 
in Ma07). 



Due to the fibres head-size and distribution, fibres cannot be 
placed closer than 25" to one another, and therefore an accurate 
and individual sky subtraction is not possible since information 
on each star's local background is not available. In addition to 
the fibres positioned on M33 stars, a set of supplementary fibres 
were added to our configurations to obtain an average estimation 
of the sky contribution. Subtracting this component allows us to 
remove both the Earth atmosphere effects and any contribution 
in the M33 field. However, the local nebular contamination will 
still be present in the final reduced stellar spectra, and must be 
considered in the analysis. The nebular contamination mainly 
affects the Balmer lines, and usually leaves the rest of relevant 
lines for morphological classification and quantitative analysis 
unaffected. 

Additionally, we were able to make use of unpublished 
spectra with a resolution of ~1.8A of M33 Var B and C and 
J013350.12+304126.6 obtained on 2003 November 30 with the 
DEIMOS spectrograph at the W. M. Keck observatory; a discus- 
sion of data acquisition and reduction may be found in Cordiner 
et al. J2008b . 

Finally, we also employed publically available optical and 
mid-IR photometric data to study the full M33 candidate LBV 
census of Ma07, supplemented with the stars in Table 2. These 
data consist of the optical photometric variability surveys of 
Shporer & Mazeh (2006) and Hartman et al. (120061) as well as 
the single epoch study of Massey et al. (2006) and the multi- 
epoch mid-IR surveys of McQuinn et al. ( 120071) and Thompson 
et al. ( 120091 ): these data are analysed in Sect. 4. 



3. Spectroscopy 

Ma07 divide their candidate LBVs into five subsets. Of these, 
the P Cygni and Ofpe/WNL stars are well defined spectral and 
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Table 1. Summary of WHT & Keck spectroscopic target list 



LGGS J# 


Aliases 


Classification 


Telescope 


V 


Last observed 


013242.26+302114.1 


- 


Iron star 


WHT 


17.44 


2006 Sept" 


013324.62+302328.4 


- 


Iron star 


WHT 


19.58 


2006 Sept" 


013350.12+304126.6 


UIT212 


Iron star 


WHT, Keck 


16.82 


1993-5* 


013406.63+304147.8 


UIT301 


Iron star 


WHT 


16.08 


200 l c 




[HS80] B416 










013357.73+301714.2* 


- 


B hypergiant 


WHT 


17.39 


2006 Sept" 


013416.44+303120.8* 


- 


B hypergiant 


WHT 


17.10 


2006 Sept" 


013422.91+304411.0* 


- 


B hypergiant 


WHT 


17.22 


2006 Sept" 


013424.78+303306.6 


- 


B hypergiant 


WHT 


16.84 


2006 Sept 


013429.64+303732.1* 


- 


B hypergiant 


WHT 


17.10 


2006 Sept" 


013339.52+304540.5 


B517 


P Cyg LBVc 


WHT 


17.50 


2006 Sept" 


013341.28+302237.2* 


[HS80] 110-A 


P Cyg LBVc 


WHT 


16.28 


2006 Sept" 


013416.07+303642.1* 




P Cyg LBVc 


WHT 


17.95 


2010 Oct'' 


013309.14+303954.5 


UIT045 


Ofpe/WNL 


WHT 


17.91 


1993-5" 


013327.26+303909.1 


UIT104 


Ofpe/WNL 


WHT 


17.95 


1993-5* 


013509.73+304157.3* 


Romano's star 


Ofpe/WNL 


WHT 


18.04 


2010 Oct'' 


013349.23+303809.1* 


VarB 


LBV 


Keck 


16.21 


2005 e 


013335.14+303600.4* 


VarC 


LBV 


WHT, Keck 


16.43 


2003-5" 


013416.10+303344.9 


UIT 341 


LBVc 


WHT 


17.12 


1993-5*' 



[HS80] B526 



Summary of of our WHT and Keck target stars, common aliases, classification following the ammended scheme of Ma07 (Sect. 3; noting 
that Romano's star should be now be considered a bona fide LBV), LGGS V band magnitude (observations made between 2000 November - 
2003 February) and date of most recent spectroscopic observations in the ~4000-5000A window ( "Ma07, ''Massey et al. 1996, Tabrika et al. 
2005, rf Neugent & Massey [20 1 1 1 g Viotti et al.|2006, ^Monteverde et al. 1996). Stars for which contemporaneous historical observations of the 
~6000-7000A window are also available in the literature are indicated by an asterix. 



evolutionary classifications, while the bona fide LBVs such as 
M33 Var B and C are also clearly identifiable; we adopt these 
categories unchanged. However the remaining groupings - the 
hot and cool LBV candidates - are more problematic since such 
nomenclature implies a specific evolutionary phase for what are 
likely heterogeneous sets of objects. 

Motivated by the appearance of M33 Var C in photometric 
minimum, Ma07 apply a classification of hot candidate LBV to 
those stars showing strong emission in both the Balmer series 
and low excitation metallic transitions such as Fen, although 
they recognise that such criteria will also include stars demon- 
strating the B[e] phenomenon (e.g. Lamers et al. 1998). This 
is intriguing, given the still uncertain physical relationship be- 
tween both classes of star and the current lack of placement of 
the latter in a coherent evolutionary scheme (see e.g. Zickgraf 
2006 for a review of this topic). Unfortunately, the limited spec- 
tral resolution and wavelength coverage of the existing spectra 
of these stars prevent us from inferring the geometry of the cir- 
cumstellar material via line profiles (e.g. Zickgraf et al. 1985) 
or H i line ratios (Lenorzer et al. 2002) and hence distinguish- 
ing between both possibilities. Moreover the Fe n emission line 
spectra of LBVs are also highly variable (e.g. AG Car, HR Car, 
R7 1 and R127; Stahl et al. 120011 Sze ifert et al. [19961 Wolf et 
al. [19811 [19881 Walborn et al. 120081 depending on the phase 
of the LBV excursion, preventing the determination of specific 
line-ratio classification criteria. 

Furthermore, such Fen-rich emission spectra are typically 
seen during and close to optical maximum, when LBVs are 
thought to be at temperature minimum, rather than in the hot 
phase supposed by Ma07 (cf. the above studies). Consequently, 
we prefer the historical term 'Iron star' for such objects (cf. 
Walborn & Fitzpatrick 2000) to simply reflect the spectral ap- 
pearance of the stars. 



Similar issues are associated with the cool LBV classifica- 
tion. Specifically, the range of temperatures experienced by an 
LBV during an excursion appears to show an explicit luminosity 
dependance (e.g. Wolf 1989). For example the high luminosity 
star R127 is observed to vary between an Ofpe to a late B-/A- 
hypergiant classification (e.g. Walborn et al. 120081 ) whereas the 
low luminosity object HD 160529 only transits between B8 to 
A9 Ia + classifications (Stahl et al. 2003); thus HD 160529 would 
be classified by Ma07 as a cool LBV in both hot and cool phases. 
Given that the examples of this classification found within M33 
are consistent with B- and early A-hypergiants (e.g. Clark et al. 
120101 120T21 Walborn & Fitzpatrick 1990), we favour a simple 
morphological classification of blue hypergiant (BHG) for these 
stars. 

Following this approach we suggest a classification of cool 
hypergiant for stars demonstrating the rich absorption+emission 
lines spectra associated with B324 and Var C in its cool phase 
(Sects. 3.5, 5.3.2 and 5.3.3). While a subset of such stars clearly 
may be LBVs undergoing excursions, some, such as B324 and 
the yellow hypergiants within Wd 1 (Clark et al. 120101 ) show 
long term stability over decadal timescales, while others, such as 
p Cas, instead transit to still cooler temperatures during outburst 
(Lobel et al. [2003 ). While we recognise that this requires the 
re-assessment of the cool LBV candidates presented by Ma07, 
as with the iron star and BHG classifications, this terminol- 
ogy relies solely on spectral morphology and carries no explicit 
evolutionary connotation, which may require additional multi- 
wavelength/epoch observations to determine. 

3. 1 . Iron stars 

We obtained spectra of four iron stars (Table 1, Figs. 1 and 
A.l). Comparison of these data to previous observations re- 
vealed that none of these stars unambiguously demonstrated be- 
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haviour characteristic of the LBV phenomenon, although signifi- 
cant spectral variability was undoubtedly present in two of them. 
This is illustrated by the spectra presented in Fig. 1, where for 
J01 3242.26+3021 14.1 we see a clear change in the emission line 
ratio of the Balmer series in the ~4 year period between observa- 
tions, with a decrease (increase) in strength of H/3 (Hy). Similar 
changes in the Balmer series over an identical time period are 
also seen in J013324. 62+302328. 4, with a dramatic decrease in 
the strength of H/3 between the two observations. 

A number of spectra of J013406.63+304147.8 are presen t 
in the literature, dating from from 1993-5 (Massey et al. 2006), 
with subsequent spectra published by Shemmer et al. ((2000), 
Sholukhova et al. d20TTb and Fabrika et al. d2005l . While sig- 
nificant multiyear gaps exist between these observations, com- 
parison of these to our new spectra suggests that there has been 
a lack of conspicuous variability over this ~17 yr period, with 
only minor changes in the line profile of Hei 4471 A and the 
anomalously strong Fe n 4556A line at a single epochQ- 

Finally, despite a >15 yr baseline, we see no gross vari- 
ation in spectral morphology between the <1995 spectrum of 
Massey et al. (fl996l l and our 2003 Keck and 2010 WHT spec- 
tra of J013350. 12+304126. 6beyond minor changes in emission 
line strength. 

We also provide the first observations of three of these stars 
in the 6000-7000A window which encompasses Ha. In each 
case it is strongly in emission - as expected given the behaviour 
of the higher Balmer transitions - with a single peaked non- 
P Cygni profile demonstrating strong emission wings. Nebular 
[Nn] and [Sn] is present in two of the four stars; given the 
comparative strength of the latter relative to the former we sus- 
pect this arises from a spatially coincident H n region rather than 
CNO processed circumstellar ejecta. Remaining emission lines 
may predominantly be attributed to He i and Fe n. 

3.2. B hypergiants 

We have made new spectroscopic observations of five BHGs 
from Ma07 (Table 1, Figs. 1 and A. 2). Historically, a classical 
temperature diagnostic for such stars is the Hei 4471A:Mgn 
4481 A line ratio (Walborn & Fitzpatrick 1 19901 1, with the 
Mgn absorption strengthening relative to Hei for later spec- 
tral types. Comparison of our spectra to those of Ma07 
suggests that this ratio has varied in two of the five stars 
examined; J013416.44+303120.8 (illustrated in Fig. 1) and 
JO 1 3422.9 1 +3044 11.0. Utilising this diagnostic, we would esti- 
mate formal spectral types of ~B8 (2006) and ~B2.5 (2010) for 
J013416.44+303120.8 (aided by comparison to examples pre- 
sented in Monteverde et al. [T996l l and ~A0 (2006) and -B5-8 
(2010) for J013422.91+304411.0, with the absence of Hei ab- 
sorption in the spectrum presented by Ma07 suggesting a later 
spectral type in 2006 than 2010 (e.g. by comparison to IFM-B 
1330; Monteverde et al. [T996] l. 

However, we caution against concluding that a secular evo- 
lution in temperature - as might be expected for an LBV ex- 
cursion - has occured on this basis alone. Firstly, we note that 
B-supergiants are prone to short period pulsational instabilities, 
as illustrated by examples found with Westerlund 1 (Clark et al. 
120101 see also Sect. 4) and we cannot rule such an explanation 

4 Foreshadowing Sect. 5.2, no characteristic LBV like excursions are 
apparent in the available photometry, although the presence of a low 
level periodic photometric modulation (Shemmer et al. 20001 and spec- 
troscopic radial velocity shifts lead Sholukhova et al. ( 20041 to propose 
that J013406.63+304147.8 is a short period binary. 



out on the basis of the current limited dataset. Moreover detailed 
quantitative non-LTE model atmosphere analysis of Galactic B- 
hypergiants by Clark et al. (2012) found that the behaviour of 
this line ratio may also depend on the location of the transition 
region between photosphere and wind; thus it does not serve as 
an unambiguous temperature diagnostic. Unfortunately, the low 
S/N ratio of the available spectra precludes us employing more 
robust diagnostics such as the Si m/Si n ratio. 

Line profile variability (LPV) is also observed in the Balmer 
transitions of these stars. This is particularly notable in the 
H/3 emission line of J013357.73+301714.2 (a weakening of the 
P Cygni profile from 2006), J013422.91+304411.0 (transition 
of P Cygni to pure emission profile), J103424.78+303306.6 
(weak P Cygni emission transitioning to absorption) and 
J013429. 64+303732.1 (the evolution of a single peaked to a P 
Cygni profile superimposed on a broad base), while variability is 
also observed in the higher transitions of J01 3357.73+3017 14.2, 
J013416.44+303120.8 and J013429.64+303732.1. Comparable 
LPV has been observed in the (wind) emission lines of Galactic 
B super-/hyp er-gia nts (e.g. Clark et al. 120101 Morel et al. [2004 
Kaufer et al. [19961 ) and LBVs (Stahl et al. 120051 . In the former 
it occurs over short (~days) timescales and has been associated 
with stochastic wind structure, while in the latter it is associ- 
ated with the LBV cycle; the current lack of sampling precludes 
us from unambiguously distinguishing between either possibil- 
ity, although the lack of other associated spectral changes argues 
against an LBV identification. 

Unfortunately, the comparatively low S/N and resolution of 
the Ha profiles presented in Ma07 precludes comparison to our 
data but we note that all five stars present an identical profile 
consisting of a narrow emission line superimposed on a broad 
base (Fig. A.2), similar to those observed in the mid- to late-B 
hypergiants in Westerlund 1 (Clark et al. 120101 ). Finally nebular 
[N n] and [S n] emission of ambiguous origin is present in four 
of the five stars observed although, as with the iron stars, we 
suspect this is not associated with circumstellar ejecta. 

3.3. P Cygni stars 

Ma07 identified four counterparts to the Galactic LBV P Cygni 
within M33, with spectra dominated by narrow emission profiles 
in the H i and He i lines; we have made observations of three of 
these (Figs. 1, 2 and A.3). As with the BHGs, LPV is present 
and observed in the H/3 line of J013341.28+302237.2 and the 
Hei lines of JO13416.O7+3O3642.10. Discovery and subsequent 
follow up spectra of both stars extend the baseline of observa- 
tions to > 17 yr for the former (Monteverde et al. 1 19961 Ma07) 
and > 18 yr for the latter (Spiller [T9921 Corral [19961 and Kehrig 
et al. 1201 II ). but provide no evidence for (LBV driven) secu- 
lar evolution during these intervals. For completeness we note 
that Urbaneja et al. (2005) prefer a classification of Bl Ia + for 
J013341.28+302237.2 (=[HS80] 1 10A) which they find to be a 
highly luminous object with a high mass loss rate in comparison 
to Galactic examples (Clark et al. 2012); we return to this point 
in Sect. 5.3.3. 

In contrast, comparison of the available dataset of 
J013339. 52+304540. 5 indicates dramatic long term variability 
(Fig. 2). Crowther et al. (119971 ) first studied this star in detail, 
finding that there was little evidence for photometric variabil- 



5 Potentially significant nebular contamination of 
J013416.07+303642.1, revealed during attempts to model this 
star, prevents us from drawing conclusions as to variability of the 
Balmer series. 
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ity between the mid-1970s to 1990s and provided spectroscopy 
from 1991, 1992 and 1995 which likewise indicated little or no 
evolution. Further low resolution spectra were subsequently ob- 
tained in 1992-4 (Corral et al. [1996b . 1993 (Massey et al.[1996) 
and 2006 (Ma07). Subject to the comparatively low resolution 
and S/N of the spectra, we find no evidence of variability be- 
tween 1991-1992, but in 1993 the Hei 4471A line appeared 
to be in absorption, before reverting to (P Cygni) emission by 
1995. However, our new spectrum indicates a dramatic change 
in morphology in the intervening period, with the disappearance 
of He i, He n and N n emission, a substantial weakening of the 
Balmer lines and the development of weak photospheric features 
in e.g. Sim. The spectrum presented by Ma07 reveals that the 
evolution had commenced prior to 2006, demonstrating an ab- 
sence of emission from high excitation species, albeit accompa- 
nied by significantly stronger emission in the Balmer lines than 
observed in either preceding or subsequent spectra. 

Model atmos phere analysis of J013339.52+304540.5 by 
Crowther et al. (fT997l > suggests T eff ~ 23 kK, M-4.9 x 
1O~ 5 M yr~' and v M ~ 275 kms -1 at this time (with the authors 
suggesting the H/He~0.8 ratio was indicative of a post LBV 
phase). We may qualitatively attribute the change in spectral 
morphology to a decrease in both stellar temperature and wind 
density, the former indicated by the disappearance of high exci- 
tation features such as He n and the latter by the weakening of the 
wind emission lines. In principle one might hope to obtain an es- 
timate of the temperature from the photospheric lines present. In 
comparison to J013341.28+302237.2 (r eff ~21 kK), the lack of 
a strong C m+O n complex at ~4650A, as well as the presence 
of Si m 4552, 68 and 75 A absorption is suggestive of a later 
classification (~B2; r e jf~18 kK), although we cannot exclude 
a somewhat earlier classification if the star is C-depleted (not- 
ing that subject to the caveats in Sect. 3.2, the absence of Mg n 
4481 A argues against a cooler temperature). Unfortunately, the 
current spectra are of insufficient S/N to be able to utilise the 
Si n and iv transitions to obtain a more accurate temperature es- 
timate; consequently we are currently obtaining suitable spec- 
tra (and contemporaneous photometry) to permit such a detailed 
model-atmosphere analysis. 

3.4. Ofpe/WNL (=WN9-11h) stars 

Ofpe/WNL (=WN9-llh) is a well defined stellar classifica- 
tion, although we note that at the low temperature extreme 
such stars share signficant morphological similarities with the 
P Cygni stars. This is well illustrated by a comparison of the 
spectra of J013416.07+303642.1 (P Cygni type; Fig. A.3) and 
Romano's star (Ofpe/WNL type; Fig. A. 4), both with one an- 
other and also with a (degraded) spectrum of P Cygni itself; 
we return to this topic in Sect. 5.3.1. Of the 12 such stars 
curently identified within M33 (Ma07, Neugent & Massev l20TlT l 
we have made observations of three; J013309. 14+304954.5, 
J013327.26+303909.1 and J013509.73+304157.3 (henceforth 
Romano's star). 

Of these, minor changes in the strength of the wind emis- 
sion lines appear present in J013327. 26+303909 (not shown 
for brevity); such variability appears common amongst Galactic 
WRs with a number of mechanisms proposed to explain this (e.g. 
binarity, pulsations, wind structure). Unfortunately, the spectrum 
of J013309. 14+303954.5 (=UIT045) is of rather poor S/N; de- 
spite this it would appear that the blend of N n+He i+He n emis- 
sion that is prominent between 4600-4700 A in the discovery 
spectrum (<1995: Massey et al.[l996) is now absent. Following 



the preceding dicussion, the spectrum now bears close resem- 
blance to the P Cygni-type stars. Nevertheless, further obser- 
vations of both both stars will be required to first validate and 
subsequently constrain the nature of the (putative) variability. 

In contrast, Romano's star is a well known photometric 
(AV ~ 16-18.6) and spectroscopic variable, demonstrating 
a ~WN8h spectrum in visual minimum and as least as late 
as WNllh in visual maximum (Maryeva & Abolmasov 2010 
Polcaro et al l20m Sholukhova et al. |20TT1FI Indeed following 
these recent studies it appears clear that it should be considered 
a bona fide LBV, a classification we adopt for the remainder of 
the paper. 

Our spectrum implies a ~WN10h classification, compared 
to a ~WN9h spectral type only 7 months earlier; consistent with 
the star simultaneously undergoing a secular photometric bright- 
ening. We highlight the close similarity to P Cygni at this time; 
the sole observational differences - the presence of weak Hen 
4686A and the lack of P Cygni profiles in the Balmer series - 
consistent with a somewhat higher temperature (Sect. 5.3.1). 

3.5. Bona fide LBVs 

Finally, we turn to the three remaining stars - J013416.10- 
303344.9 and M33 Var B and C - which Ma07 denote as LBVs 
(inevitably without spectroscopic classification criteria due to 
their variability; Table 1 and Fig. A.5). The 4000-5200A spec- 
tra of J013416.10+303344.9 and M33 Var C obtained in 2010 
and the 2003 spectrum of M33 Var B share very similar mor- 
phologies, being dominated by P Cygni Balmer and Fe n emis- 
sion lines, which are strongest in the latter star. P Cygni Hei 
emission lines are present in M33 Var B and C, but are in ab- 
sorption in J013416. 10+303344.9. The spectral morphology of 
M33 Var B and C, when observed in 2003 and 2010 (respec- 
tively), appear to resemble composites of the P Cygni and iron 
star classifications (Fig. A.5), suggesting that the physical prop- 
erties of these stars at the times of these observations were inter- 
mediate between these subtypes. Indeed, the spectrum of M33 
Var B is strongly reminiscent of P Cygni (Fig. A.5) - although 
the presence of Fe n emission implies a cooler, denser wind than 
that star. 

In contrast, the 2003 spectrum of M33 Var C is of a very 
different character, being dominated by a wealth of low excita- 
tion metallic absorption lines, with only H/3 weakly in emission. 
Comparison to the 2010 observation and the highly luminous 
F0-5 Ia + hypergiant B324 (Fig. 2; Monteverde et al. 1996) indi- 
cates a significantly later spectral type at this epoch. While the 
Ha line is saturated in the 2003 spectrum (as it is in the other 
2 DEIMOS/Keck spectra), weak P Cygni emission is present in 
the 40 and 70 multiplets of Fe n; again similar to B324 (Sect. 
5.3.3). 

Regarding variability, the candidate LBV 
J0134 16.10+ 303344.9 was ob served by both Monteverde 
et al. j 19961 1 and Massey et al. (1996) between 1993-5; com- 
parison of these data to our new spectrum show a broadly 
similar morphology, although the Balmer lines appear to have 



6 Limited low resolution spectroscopy from 1992, coincident with the 
absolute visual maximum (Szeifert 1996 ) is suggestive of a somewhat 
later B-supergiant classification. These appear not to be dominated by 
the Fe n emission that characterises e.g. AG Car, HR Car or R 127 during 
visual maximum (Stahl et al. 2001 Szeifert et al. 2003 and Walborn et 
al. 2008 ), with the apparent lack of low excitation metallic photospheric 
lines suggesting an earlier classification than the F-supergiant spectrum 
demonstrated by e.g. M33 Var B (Ma07). 
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strengthened in this period. In contrast, the LBV classifications 
of M33 Var B and C are long established (e.g. Szeifert et al. 
119961 Viotti et al. 120061 and refs. therein). Ma07 present the 
most recent high resolution blue end spectra of M33 Var B 
that we are aware of (from 1993 and 2005); our spectrum 
differs notably from both of these, although in the absence of 
a contemporaneous lightcurve or better temporal spectroscopic 
coverage we do not discuss this further. Finally, while clearly 
highly variable, we defer analysis of the spectral evolution of 
M33 Var C to Sect. 5.3.2, where it is discussed in the context of 
the long term photometric behaviour of the star. 

4. Photometry 

Photometric variability is a defining property of LBVs, and in 
addition to our new spectra, we may utilise existing photomet- 
ric surveys (e.g. Hartman et alj2006 Shporer & Mazeh 2006, 
McQuinn et al. [20071 Thompson et al. [2D09l > of M33 to further 
constrain the nature of the candidate stars (Ma07, supplemented 
with objects listed in Table 2). 

4. 1 . Optical photomety 

Ma07 briefly reviewed the optical photometric data available 
for the candidate LBVs within M33, identifying variables both 
by comparison of historical to Local Group Galaxies Survey 
(LGGS; Massey et al. 2006) observations and also via the multi- 
epoch photometric surveys of Hartman et al. (2006) and Shporer 
& Mazeh (2006). Expanding on this, we retrieved and exam- 
ined the individual lightcurves of the 15 stars flagged as vari- 
ables by these studies. Of these, only one star demonstrated pre- 
viously unidentified, coherent long term variability - the iron star 
J013459.47+303701.0 (Fig. 3). While this is potentially consis- 
tent with an LBV identification we note that the star became 
bluer as it brightened, whereas the opposite behaviour might be 
expected, with LBVs typically cooling as they brighten. 

Of the known LBVs, the 2000-2003 lightcurve of M33 Var C 
(Shporer & Mazeh 2006) is of particular interest (Fig. 4), since 
it indicates a ~0.8 mag. brightening between 2001-2 July, fol- 
lowed by an apparent plateau with V ~ 15.4 mag (comparable 
to the maximum light of the two previous maxima; Humphreys 
et al. 1 1 9881 ). We return to discuss the combined photometric and 
spectroscopic dataset for this star in Sect. 5.3.2. 

After excluding known (photometrically variable) LBVs 
such as M33 Var C from the sample, apparently aperiodic 
photometric variability of <0.4mag. amplitude was associated 
with the remaining objects (e.g. J013426.1 1+303424.7; Fig. 2fl 
Likely driven by pulsations, quasi-periodic photometric variabil- 
ity over timescales of tens to hundreds of days appears ubiqui- 
tous amongst supergiants of spectral type B and latter - the so 
called a Cygni variables (van Leeuwen et al. 119981 Clark et al. 
120101 ) - as well as LBVs (Lamers et al. 119981 ) and it is tempting 
to attribute the variability to this cause (with a combination of 
changes in the pulsational period and poor temporal sampling 
accounting for the lack of identifiable periodicity). Finally, we 
note that despite rather large amplitude (AV ~ 0.8mag.) his- 
torical photometric and spectroscopic variability (Ma07; Fig 1), 

7 The iron stars J013235.25+303017.6, J013242.26+302114.1, 
J013300.02+303332.4, J013406.63+304147.8 and 

J013426.1 1+303424.7, the Ofpe/WN9 stars J013309.14+304954.4, 
J013327.26+303909.1 and J013353.60+303851.6, the B hypergiants 
J013355.96+304530.6 and J013416.44+303120.8 and the P Cygni 
candidate J013416.07+303642.1. 



the iron star J013242. 26+3021 14. 1 was not flagged as variable 
by either survey, nor were a number of the spectroscopic vari- 
ables identified in this worlfl including the strong LBV candi- 
date J013339.52+304540.5 (Sect 3.3 and Fig. 2; although pre- 
empting the following section, the BHG J013429.64+303732.1 
is a mid-IR variable). 

We may also utilise the LGGS photometry to construct 
colour manitude plots for the candidates of Ma07, which we 
show in Fig. 5. We can clearly identify the blue plume of 
highly luminous blue supergiants that Massey et al (2003 de- 
scribe. However, plotting both (B - V) and (V - I) colour in- 
dices demonstrates that a number of individual objects such 
J013340. 60+304137.1 show rather discrepant colours, unfortu- 
nately making dereddening individual stars problematic. 

Nevertheless, of the objects plotted, the brightest, 
J013350.92+303936.9 (=UIT 218, which Ma07 label as a 
hot LBV candidate) is in fact the spatially extended nucleus of 
M33, which Massey et al. ([1996) show has an F supergiant like 
spectrum (and hence we do not consider it further). The next 
brightest, J013355.96+304530.6 (=B324) appears to be both 
intrinsically bright and possesses the latest spectral type of any 
object considered here - we return to this in Sects. 4.2 and 5.3.3. 

Of the remaining objects from Ma07, the Ofpe/WNL stars 
appear to be the faintest subtype, with V > 17.5, although there 
is some overlap with both the P Cygni and B hypergiants. With 
V < 18 and < 17.5 respectively, these subtypes appear system- 
atically brighter (as might be expected from the larger visual 
bolometric corrections anticipated for the Ofpe/WNLs). Given 
their intrinsic variability, bona fide LBVs span a wide range of 
magnitudes, with Romano's star having been found to be at least 
0.6 mag fainter in the V band than shown here (e.g. Sholukhova 
et al. 1201 II ). Finally, and somewhat unexpectedly, the iron stars 
are found to span the largest magnitude range (AV ~3.5 mag); 
we return to the nature of these objects in Sect. 5.2. 

4.2. Mid-IR photometry 

M33 has been the subject of extensive mid-IR studies, with six 
epochs of Spitzer/IRAC (3.6-8.0/vm) observations between 2004 
January and 2006 February (unfortunately not contemporaneous 
with the LGGS photometry or the spectroscopy of Ma07 and 
this work). Point source catalogues derived from these data have 
been compiled by McQuinn et al. (120071 ) and Thompson et al. 
(2009) enabling identification with the source list of Ma07; a 
task undertaken by the latter authors. To this list we were able 
to add the additional candidates listed in Table 2 and cross cor- 
related the expanded census against both catalogues (as well as 
adopting fluxes for four sources kindly provided by Valeev 20 1 1 , 
priv. comm.). 

We found a total of 29 sources had a detection at one or more 
wavebands. Of these 18 were from the catalogue of Thompson 
et al. ( 120091 including three BHGs previously unconsidered by 
these authors), four supplied by Valeev et al. (2009) and seven 
new counterparts unique to McQuinn et al. d2007l l. Of the 16 
stars from the catalogue of Ma07 identified as having mid-IR 
counterparts by Thompson et al. (2009), 13 were found in the 
catalogue of McQuinn et al. (2007). These showed systematic 
offsets of 1113.5 ~0.15 mag and 1114 5 ~0.2 mag in the sense that 
the magnitudes in Thompson et al. (2009 ) were brighter; a corre- 
sponding correction has not been applied to the values presented 



8 J013324.62+302328.4, J013341.28+302237.2, 
J013422.91+304411.0, J013357.73+301714.2, 
J0134234.78+303306.6 and J013429.64+303732.1 
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Fig. 4. V band photometry of M33 Var C from Shporer & Mazeh 
(2006), with the timing of the spectroscopic observation dis- 
cussed in Sect. 3.5 indicated. 

show significantly greater offsets than described her^l in such 
cases we adopt the relevant fluxes from Thompson et al. (2009 ). 

Broken down by nature of sourc43 14 of 18 iron stars, 9 of 
the 13 BHGs, one each of the 1 1 Ofpe/WNL and four P Cygni 
candidates, all four LBVs and the sole cool hypergiant, B324, are 
detected in one or more wavebands. Fluxes for those sources not 
listed by Thomspon et al. (2009; their Table 3) are given in Table 
3 of this paper. Only three stars are found to be variable; M33 Var 
C, J01 3429.64+303732.1 (Thompson et al. [20091 and Romano's 
star. The lightcurve of M33 Var C (Fig. A. 3 of Thompson et al. 
2009) is particularly interesting, since fortuitously the Spitzer 
observations commenced at the point when the optical observa- 
tions of Shporer & Mazeh (2006; Fig. 4) ended and indicate a 



in Table 3 but has been applied in Figs. 6 and 7 to enable direct 
comparison. No obvious offset was found with respect to 8jt/m 
fluxes, although we caution that individual sources in all 3 bands 



9 J013333.22+303343.4 at 8pm, Var B at 4.5pm, 
J013350.92+303936.9 at 3.6pm and J013416. 10+303344.9 at 4.5pm 

10 Excluding the nuclear cluster source J013350. 12+304126.6, reclas- 
sifying Romano's star as a bona fide LBV and adopting a classification 
of iron star for N45 901, N93351 andN125093 (Table 2). 
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~0.8mag fade over a period of ~400 days commencing early in 
2005; we return to this in Sect. 5.3.2. 

A number of different physical processes might be expected 
to contribute to the ~3-8//m fluxes of the sources, including pho- 
tospheric emission (potentially from a putative cool companion), 
free-free continuum emission from a dense wind and thermal 
emission from a dusty circumstellar environment (either result- 
ing from dust condensation in a wind shock in a colliding wind 
binary or in red supergiant/LBV ejecta). Given this, the current 
uncertainty in the parameters of the underlying stars (and hence 



the stellar contribution to the composite spectral energy distribu- 
tion) and the potential difficulty in accounting for the contribu- 
tion of line and molecular (e.g. PAH) emission in the continuum 
bands, quantitative analysis of individual sources appears pre- 
mature. Nevertheless, we may attempt to determine the likely 
nature of the mid-IR emission from different objects in a statis- 
tical sense utilising various different mid-IR colour/magnitude 
and colour/colour diagnostic plots that have been developed for 
such evolved massive stars (e.g. Hadlield et al. 2007, Bonanos et 
al. |2009l . 



10 



J. S. Clark et al.: On the nature of candidate luminous blue variables in M33 



6 



10 

CD 

CO 



A 



,2 







Object X m 



J01 3426.1 1+303424.7 • 

,mt M33 Var A 



• LMC sgB[e] 

• LMC LBV+WN9-llh 

. LMC WN 

(excl. WN9-llh) 

• LMC WC 



imqqoa co onoooo a ° LMC 0BA superfliant 

J01 3324.62+302328,4 a . JQ1 3459.47+303701 .9 

J01 3500.30+3041 50.9* " 



J01 3333.22+303343.4 . 
J01 3350.1 2+3041 26.6 

V. 

V: .-. 



■ N93351 
,N1 25093 



.J01 3442.1 4+30321 6.0 



Romano's star 



J01 3248.26+303950.4 



J01 3424.78+303306.6 




J01 341 0.93+303437.6 (Var 83) 
^013242.26+302114.1 J01 3349.23+303809.1 (Var B) 

U01 3341 .28+302237.2 ■ 
013335.14+303600.4 (Var C) 



G 



2 



[4.5] [8.0] 
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Following Thompson et al. we initially utilise the [3.6]-[4.5] 
vrs. M4 5 colour/absolute magnitude plot (Fig. 6 - employing a 
distance modulus, ^=24.92; Bonanos et al. I2006I I. The stellar 
population of M33 is predominantly located in the lower left 
of the plot, with a plume extending to larger M4 5, that likely 
represents a combination of foreground and intrinsically lumi- 
nous (composite) sources. A population of dusty, (extreme) AGB 
stars appear to dominate the region of the diagram delineated 
by M45 > -11 and approximately [3.6]-[4.5]>0.5 (Thompson 
et al. 120091) . Finally, a population of massive, evolved stars are 
located in the upper right corner. Following these authors and 



Khan et al. d201 II) . we plotted selected Galactic objects here for 
comparison, noting that the [3.6]-[4.5] colour indices of these 
stars suggest the presence of hot dust (with a simple black body 
temperature, r BB ~ 350 - 1000 K). 

We may supplement this with the [ 3.6]- [4.5] vrs. [4.5] -[8.0] 
colour colour plot of Bonanos et al. (f2009; Fig. 7). The nature 
of the emission from sources in this diagram is best demon- 
strated by plotting the massive stellar population of the Large 
Magellanic Cloud (excluding the cool supergiants; Bonanos et 
al.). The OBA super-/hypergiant and WRs form a continuous se- 
quence of increasing mid-IR colour excess, with the latter redder 
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than the former. This may simply be understood as the increas- 
ing contribution from free-free emission, which on average is 
more pronounced for the WRs due to their denser winds. LBVs 
and Ofpe/WNL stars are also found along this sequence, with 
an additional contribution from circumstellar dust explaining the 
outliers. Finally, stars exhibiting the B[e] phenomenon form a 
distinct group, with their very red colours the result of the hot 
dust that characterises such objects. 

In the colour magnitude plot (Fig. 6), the majority (18) of 
the stars with mid-IR detections appear to be broadly co-located 
with the underlying stellar population of M33 (approximately 
bounded by M 45 > -11 and [3.6]-[4.5]<0.5); these consist of 
representatives of each of the ammended morphological classi- 
fications of Ma07 (bar the single cool hypergiant, B324). In the 
absence of quantitative modeling, the nature of the IR emission 
of these stars is uncertain, although we suspect it is due to a com- 
bination of wind (free-free) and photospheric emission. Indeed 
the three mid-IR variables - M33 Var C, Romano's star and the 
BHG JO 13429.64+3032 16.0- are located in this region and none 
show correlated mid-IR colour/magnitude changes (e.g. in the 
sense that they become redder as they brighten), as might be ex- 
pected if the emission originated in a variable dusty circumstellar 
environment. 

Unfortunately, only seven of these stars have 8/im detec- 
tions (Fig. 7). Of these, Romano's star and M33 Var C are co- 
located along the OBA supergiant/WR sequence; suggestive of 
wind+photospheric emission (although Valeev et al. 120091 at- 
tributes the IR emission to hot dust, highlighting the difficulty in 
interpretation). The remaining objecto are all displaced from 
this locus by virtue of their comparatively strong 8^m emission, 
implying the presence of warm circumstellar dust. Note how- 
ever, that the absence of a pronounced 8^/m flux does not mean 
that the remaining stars lack dusty ejecta, since emission from 
LBV nebulae typically peaks at significantly longer wavelengths 
(e.g. Egan et al. 120021 Clark et al. 120031 Voors et al. 120001 . 

A further 11 stars are found outside this region of the 
colour magnitude plot. Of these, the mid-IR colours of the iron 
star J013242. 26+3021 14.1 are suggestive of wind+photospheric 
emission (Fig. 7) and while no 8/mi detection is available for 
J013355.96+304530.6 (=B324), we suspect a similar origin for 
its 3.6 - 4.5/mi flux, given both its intrinsic extreme luminosity 
and strong wind emission lines (Sect. 5.3.3). 

This then leaves nine mid-IR luminous stars (M4 5 < -11), 
all of which are iron starQ. These are found in a sequence 
of increasing M45 and [3.6]-[4.5] colour (Fig. 6) that termi- 
nates amongst highly luminous evolved stars with rich dusty 
circumstellar environments such as the Galactic RSGs NML 
Cyg (Schuster et al. 120091 an d VY CMa (Smith et al. I200T1) . 
the LBV 77 Ca r (Sm ith et al. [20021 . the YHG IRC +10 420 
(Blocker et al. 119991 1 an d, with in M33, the cool transie nt M3 3 
Var A (Humphreys et al. 120061 ) and Object X (Khan et al EOTlT ). 
Of these, four - J013324.62+302328.4, J013426.1 1+303424.7, 
J013459.47+303701.9 and J013500.30+304150.9 - are co- 
located with the sgB[e] stars (Fig. 7), while the colours of a 
further two - J013333.22+303343.4 and J013350. 12+304126.6 
- are also suggestive of (a smaller quantity of) hot dust. The fi- 
nal three stars - N93351, N125093 and J013442. 14+303216.0- 
also appear associated with hot dust, although their larger 8/mi 



11 The LBVs M33 Var 83 and Var B, the iron star 
J013248.26+303950.4, the BHG J013424.78+303306.6 and the P 
Cygni star J013341.28+302237.2 

12 Note that identification of the most luminous of these sources as 
M33 Var C in Thompson et al. (2009) is in error. 



fluxes imply greater quantities of warm dust is also present with 
respect to the former objects. Unfortunately, due to the lack of 
longer wavelength data we refrain from attempting to determine 
temperature(s) for the dust component(sQ although Bonanos 
et al. ( 120091 ) report dust temperatures of ~600K for those LMC 
sgB[e] stars for which 24/vm fluxes are available. 

These stars are found to span an unexpectedly wide 
range of ~3 magnitudes in the V band (Fig. 5), includ- 
ing objects that appear to be either intrinsically faint (e.g. 
J013324. 62+302328. 4) and/or subject to significant reddening 
(e.g. J013459.47+303701.9). Similar behaviour is also found 
at mid-IR wavelengths, with the nine stars spanning >2.3 
mags in the 4.5/vm bandpass. The two brightest objects - 
J013350.12+304126.6 and N125093 - are comparable to the 
cool transient M33 Var A and the Galactic YHG IRC +10 420 
(Fig. 6). However their mid-IR colours suggest differences in 
dust temperature in comparison to these objects, particularly for 
J013350. 12+304126.6 which, by virtue of its 8/vm flux, appears 
to lack a warm dust component. 

Turning to the fainter extreme and we find that the four 
faintest mid-IR sources from this subset are those that have mid- 
IR colours consistent with sgB[e] stars. With V > 18 mag. 
these are also amongst the faintest stars at optical wave- 
lengths (Fig. 5). Bonanos et al. (2009) report that the LMC 
sgB[e] stars experience greater reddening than normal OB stars; 
given the reddened optical colours of J013426.1 1+303424.7 and 
J013459. 47+303701. 9 we suspect that this may also be the case 
for these stars. An (extreme) example of this behaviour would 
then be the evolved massive star 'Object X', where extensive re- 
cent mass loss has rendered it the brightest stellar object within 
M33 in the mid-IR, but undetectable at wavelengths shortwards 
of the V band (Khan et al. 120111 ). Applying magnitude (M4.5 < 
-11.5) and colour (0.5<[3.6]-[4.5]<1.5) cuts designed to iden- 
tify massive dusty objects, Thompson et al. (120091 ) identified a 
further 16 optically faint and 9 optically undetectable sources 
within M33, and it is tempting to attribute our subset of opti- 
cally faint iron stars associated with hot dust to a low extinction 
tail of this population; we return to the nature of these objects in 
Sect. 5.2. 



5. Discussion 

A central goal of this work is the determination of the physical 
and evolutionary nature of the stellar sample. It is easy to place 
both the Ofpe/WNL and P Cygni stars at the high temperature 
extreme of an LBV excursion, with the former appearing to be 
both hotter and to have larger wind terminal velocities (for a 
comparable mass loss rate, leading to a lower wind density) than 
the latter (e.g. Clark et al. 1201 lal and in prep.). However, how 
may the remaining stars be accommodated? 

5.1. The BHGs 

Clark et al. (120121 ) studied the Galactic population of BHGs and 
found that they appeared to result from two distinct channels, 
with high luminosity examples such as Sco found in a pre- 
LBV phase and lower luminosity objects such as HD 160529 



13 Indeed, even with longer wavelength data, the study of Object X by 
Khan et al. (201 1 ) illustrates the difficulty in quantitative analysis, pre- 
senting two acceptable fits to the observed spectral energy distribution 
arising from the adoption of very different temperatures and chemical 
properties of the circumstellar dust. 
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Table 3. Mid-IR fluxes for new candidate LB Vs and BHGs 



Name 


t^lass. 


L3.6J 


L4.5J 


[5.8J 


[8.0J 


UIT005 


BHG 


16.46 


16.35 


~ 


- 


[HS80] B133 


BHG 


16.33 


16.20 






1FM-B 1330 


BHG 


15.96 


15.95 






M33 Var B 


LBV 


15.65 


15.53 


14.82 


13.20 


N93351 


iron star 


13.30 


12.64 


11.67 


10.60 


N 125093 


iron star 


12.36 


11.73 


11.05 


10.00 


Romano's star 


T T}\7 

Liiv 


14. yj 


1 A H£. 
14. /O 


1 A 70 

14. 1 L 


1 A 1 A 

14. /4 


013235.25+303017.6 


iron star 


16.32 


15.87 






013242.26+302114.1 


iron star 


13.70 


13.61 




13.18 


013309.14+302328.4 


Ofpe/WN9 


17.32 


17.01 






013350.92+303936.9 


iron star 


10.99 


11.08 




11.40 


013355.96+304530.6 (=B324) 


YHG? 


12.74 


12.52 






013357.73+301714.2 


BHG 


16.74 


16.53 






013416.44+303120.8 


BHG 


16.61 


16.66 







Top panel: fluxes from the catalogue of Thompson et al. (2009), middle panel: fluxes from Valeev et al. (2009. priv. comm. 2011) and bottom 
panel: fluxes from the McQuinn et al. ([2007 ). Note that the values from McQuinn et al. (12007 ) have not been subject to correction due to 
a systematic offset with those of Thompson et al. (2009), while the former do not provide [5.8] magnitudes. Classifications given follow the 
ammended morphological scheme of Ma07, noting that Massey et al. ( 1996 ) associate 013350.92+303936.9 (=UIT218) with the galactic nucleus; 
hence it is likely that it is a composite source. 



being post-RSG objects potentially encountering the LBV phe- 
nomenon as they evolve to higher temperatures. Given the con- 
cerns regarding both dereddening individual stars and the use 
of the He i 4471/Mgn 4481 A ratio as a temperature diagnos- 
tic (Sect. 3.2) we refrain from calculating bolometric luminosi- 
ties for individual BHGs discussed in this study. Nevertheless 
assuming negligible reddening, a bolometric correction appro- 
priate for an AO spectral classification (BCA0=-0.15mag and 
for reference BCB2=-l-59mag; Clark et al. 20053) an d a lower 
limit of V~17.5 for the range of optical magnitudes spanned 
by these stars, results in a minimum bolometric luminosity of 
log(L/L o )~5.0; consistent with both low and high luminos- 
ity/mass channels. Thus their individual placement in an evolu- 
tionary scheme will have to await systematic model atmosphere 
analysis, although of the three BHGs subject to such a study 
(Urbaneja et al. [20031 l2QTTb [HS80] 1 10-A is found to be unex- 
pectedly luminous, (Sect. 5.3.3), while the stellar+wind proper- 
ties of both UIT 005 and OB 10- 10 are similar to those of pre- 
LBV Galactic BHGs (Clark et al. [2012b . 

5.2. The iron stars 
5.2.1. LBV candidates 

While resulting in a heterogenous group of stars (Sect. 3), the 
observational criteria employed by Ma07 to define the class of 
iron stars (Sect. 3.1) are clearly successful in identifying LBV 
candidates. The candidatures ofN45901, N 93351 and N125093 
are discussed by Valeev et al. (2009, 2010), and are supported by 
pronounced photometric variability (N9335 1 and potentially the 
remaining stars), mid-IR excesses (N93351 & N125903; Sect. 
4.2) and the presence of [Can] AA 7291, 7323 emission (N93351 
and N 125093), which appears associated with variable/eruptive 
stars such as M33 Var A and IRC +10 420. The impressive sim- 
ilarity of the spectrum of J013350.12+304126.6 to that of the 
LBV R127, together with its pronounced mid-IR excess, flags 
it as a strong candidate even in the absence of current variabil- 
ity, while the mid-IR properties of J013442. 14+303216.0 also 
points to the presence of a cool dusty circumstellar environment. 
J013459. 47+303701. 9 also appears dusty and demonstrates co- 



herent photometric variability (Fig. 3); although it is optically 
fainter than the preceding examples this may in part be due to 
greater (intrinsic) reddening. Finally Ma07 report striking spec- 
troscopic variability in J013332.64+304127.2 (from a WNL to 
iron star morphology), although unlike the preceding examples 
this star does not appear to support an IR excess. 

Therefore, under the assumption that a subset of both the 
iron stars and BHGs are indeed bona fide LBVs, consideration 
of confirmed high and low luminosity LBVs (e.g. AG Car, HR 
Car, R71, R127 and HD 160529) allows us to infer a broad 
physical progression from high to low temperatures in terms of 
the ammended classification scheme of Ma07: 

Ofpe/WNL — > P Cygni — > iron star/BHG — > cool hypergiant 

This sequence is subject to three caveats; that (i) progress 
between subtypes is continuous (cf. the apparent hybrid spectra 
of M33 Var B; Sect. 3.5), (ii) the precise appearance of the spec- 
tra are also subject to the properties of the stellar wind, which 
are also found to vary through an LBV excursion (e.g. Groh et 
al 2009) and (iii) depending on stellar luminosity, the full range 
of subtypes may not be experienced, with lower luminosity ex- 
amples not reaching the hottest Ofpe/WNL and P Cygni phases 
(cf. Wolf [1989l . 

5.2.2. B[e] star candidates 

As discussed in Sect. 4.2 we have identified a subset of objects 
with IR excesses consistent with the hot dust that characterises 
sgB[e] stars and which are amongst the optically faintest of our 
sample (see Thompson et al. 2009 for individual SEDs). While 
this in part may be the result of a high degree of (circumstellar?) 
reddening towards certain objects, both J013500.30+304150.9 
and J013324. 62+302328.4 appear to suffer little extinction and 
hence are likely to be intrinsically faint (Fig. 5). Given the ab- 
sence of photospheric features, determining bolometric correc- 
tions for both these stars is difficult. However, adopting values 
from Clark et al. d2005bl) appropriate for the range of spectral 
types suggested for Magellanic Cloud sgB[e] stars (B0 to B8; 
Zickgraf et al. 119861) we find, prior to applying an uncertain 
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reddening correctiorQ log(L„/L ) -4.1-5.1 and -4.2-5.2 for, 
respectively, J013324.62+302328.4 and J013500.30+304150.9. 
Such luminosities are below those of the sgB[e] stars analysed 
by Zickgraf et al. ( 1986 ) and known Galactic LBVs (Clark et al. 
2005a), but they are consistent with a lower luminosity popula- 
tion of stars demonstrating the B[e] phenomenon identified by 
Gummersbach et al. (I1995I ). 

What might be the physical nature of the B[e] star can- 
didates? As described by Lamers et al. (1998), the B[e] phe- 
nomenon occurs in both single and binary stars with a range of 
masses and evolutionary states. The lack of apparent molecu- 
lar bands in the spectra of these stars would appear to exclude 
symbiotic systems, while the mid-IR colours appear inconsis- 
tent with those of (massive) young stellar objects (e.g. Bolatto et 
al. 2007). If a subset of stars are confirmed to be of low luminos- 
ity (log(L*/L ) < 4.7), they would have corresponding progen- 
itor masses of < 15M (under the assumption they have left the 
main sequence; Meynet & Maeder 2000). While they may not 
be directly identified with the extreme AGB stars identified by 
Thompson et al. (120091 ) - since their emission spectra imply they 
must be substantially hotter - such luminosities are comparable 
to those of the dust enshrouded progenitors of SN 2008S and 
the 2008 transient in NGC300, which are found to be optically 
obscured, mid-IR bright and non variable (Prieto 2008 Prieto et 
al. 12008 Thompson et al. 120091 ). Consequently it is tempting to 
consider the possibility that the 2 classes of object may be re- 
lated, but subject to differing intrinsic absorptions; a suggestion 
first advanced for low luminosity sgB[e] stars in the LMC by 
Bonanos et al. (120091 . 

More luminous candidates would be comparable to the bona 
fide sgB[e] stars in the Magellanic Clouds. While the physical 
mechanism that leads to the production of sgB[e] stars is cur- 
rently uncertain, by analogy to classical Be stars it is suspected 
that high rotation rates are necessary for the formation of a cir- 
cumstellar disc. Both binary interactions and mergers have been 
invoked to spin up the star, while a variant of this scenario - the 
presence of a circumbinary disc produced by short lived binary 
drive mass transfer - has been invoked to explain the sgB[e] star 
Wdl-9 (Ritchie et al. in prep). Given this, it is particularly in- 
teresting that radial velocity shifts led Sholukhova et al. (2004) 
to conclude that J013406. 63+304147.8 is a short period binary, 
with mass loss in the equatorial plane leading to a sgB[e] classi- 
fication. 
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Fig. 8. Comparison of the synthetic spectrum (dashed) to the 
2010 Sept. observations (solid) of Romano's star. 



5.3. Individual objects 
5.3.1. Romano's Star 

As highlighted in Sect. 3.4, observational studies clearly reveal 
the LBV nature of Romano's star. Recent quantitative modeling 
of the star during the optical minimum of 2007/8 and the lo- 
cal optical maximum in 2005 by Maryeva & Abolmasov (2010) 
reveal significant differences between the two epochs in the 
sense that the star was both hotter and more compact during 
optical minimum, as expected for canonical LBV excursions. 
Interestingly, however, these authors report a significant differ- 
nce in bolometric luminosity between the two phases, whereby 
the star was more luminous during optical minimum; a point we 
return to later. 

Given that a central goal of our program is to quantita- 
tively follow the evolution of the stellar parameters of LBVs 

14 Ma07 give an average E(B — V) = 0.12 for blue stars within M33 
which corresponds to ~0. 15dex in log(L„); dwarfed by the uncertainty 
in the temperature dependant bolometric correction. 



through their excursions/outbursts, Romano's star provided an 
ideal proof-of-concept; moreso since our spectroscopy and quasi 
simultaneous photometry (Sholukhova et al. 1201 II) indicates that 
it appeared to be in a state intermediate between those de- 
scribed by Maryeva & Abolmasov (1201 II) . To accomplish this 
we utilised the non-LTE model atmosphere code CMFGEN 
(Hillier & Miller [19981 [19991 ) and followed the methodology 
previously employed for the spectroscopically similar stars P 
Cygni (NaiaiTo l2U0Tl ) and NGC 300 B16 (Bresolin et al. 2002), 
with the observational constraints provided by our September 
2010 blue and red spectra and the photometric values from 
Sholukhova et al. d20TTl) . 

The latter were obtained within one month of our spectro- 
spic observations and according to Sholukhova et al. (201 1), the 
object's V magnitude varied between 17.75 and 17.85, with B 
ranging from 17.65 and 17.75, resulting in an average (B — V)~ 
-0.10 (Fig 2.b of Sholukhova et al. l20TTl) . During this epoch, the 
object was moving from the hotter, minimum brightness to the 
cooler, maximum brightness phase, being closer to the former. 
Fortunately, though rather weak, the Heii/14686 line could be 



14 
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used to constrain the ionization equilibrium and therefore derive 
the effective temperature of the staiFT. 

Nevertheless, when combined with the ionization status dis- 
played by the star, the moderate S/N of our spectra hamper an ac- 
curate determination of all stellar properties. Reliable estimates 
could be obtained for the effective temperature, , stellar lumi- 
nosity, L», the wind terminal velocity v,x,, the H/He ratio and the 
mass-loss rate, M. On the other hand, stellar parameters such as 
the exponent of the velocity field (J3), degree of wind clumping 
(/v) and nitrogen and silicon abundances are subject to larger 
uncertainties, while only upper limits may be returned for ad- 
ditional elements (e.g. Fe, C, Mg). In order to investigate this 
behaviour we ran a suite of models at solar, 0.7 x solar and 0.4x 
solar metallicity and, apart from the Si lines, found no signifi- 
cant differences in either the quality of the final fit or the basic 
stellar properties, noting that higher S/N observations at maxi- 
mum brightness (cooler phase) would unambiguously constrain 
the metallicity of the object. Because of this result metallicity 
had to be assumed for our analysis. To accomplish this we uti- 
lized [O/H] estimates from the literature based on the radial po- 
sition of the star within M33. According to Romano (1 19781 1 the 
star is placed at a projected distance of 17' (4.8kpc) from the 
center of the galaxy, for which Neugent & Massey ( 201 11 1 and 
Bresolin (120021 derive [O/H] -8.25 - 8.30 and [O/H] -8.30 re- 
spectively, corresponding to 0.4xsolar. We adopted this value 
for the abundances of Na, Mg, Al, S, Ca, Ti, Cr, Mn, Fe, Co and 
Ni and assumed a distance of 964kpc to M33 (Sect. 4.20 

The results of this analysis are summarised in Table [4] while 
we compare our best-fit synthetic spectrum to observations in 
Fig. 8. As the emission strength of He H/14686 is extremely sen- 
sitive to the effective temperature in this parameter domain we 
derive an uncertainty of ~500K for the r e ff of the object. Taking 
into account the uncertainties in the distance, photometry and 
reddening we estimate an error of 0.15dex for the stellar lumi- 
nosity. Wind properties are fairly well constrained. The terminal 
velocity, Voo= 265 + 50km s , is set by the P Cygni profiles 
of the He i lines which, together with the shape of the emission 
components of Ha and HB, help constrain B. We adopted a two 
/Haw with a flat B — 6 component in the inner parts and a steeper 
B = 0.9 law in the outer wind. Despite the medium S/N of our 
spectra, the strong electron scattering wings of Ha point towards 
the presence of moderate clumping, /v=0.25. 

Given the moderate spectral resolution of our observations 
we were not able to unambiguously identify and subtract po- 
tential nebular contamination of the Hi lines introduced by 
the fibre size (1.6" diameter; Sect. 2). To compensate for this 
we placed larger weight on reproducing the Hei and higher 
Hi Balmer lines than on Ha and HB when determining M. 
Thus our current favoured model slightly underestimates emis- 
sion in the latter transitions and we estimate a likely error of 
<25% on this parameter in light of this uncertainty. We obtain 
H/He=1.5 by number with H/He=2.0 and H/He=1.2 as upper 



We highlight that over the wavelength range discussed and at the 
S/N of our observations, the Hen 4686A line is the sole viable tem- 
perature diagnostic that permits spectroscopically similar objects such 
as P Cygni (T cff ~ 18kK; Naiarro I200TV J013416.07+303642.1 and 
Romano's star (T eff ~ 26kK; Figs. A.3 and A.4) to be distinguished. 
Spectra intended for quantitative modeling of such objects must there- 
fore be of sufficient resolution and S/N to resolve this feature, while 
higher S/N observations would also permit N n/N m and Si m/Si iv ioni- 
sation equilibria to be employed. 

16 Maryeva & Abolmasov J2011b adopt a slightly smaller distance 
(847 kpc) implying /j=24.64. Correcting the stellar parameters for this 
lower distance does not materially affect our conclusions. 



and lower limits respectively. The silicon abundance, with an 
uncertainty of 0.2dex, is basically constrained by the strength of 
the Si ni/14550 - 4575 lines. Our derived nitrogen abundance re- 
produces satisfactorily the main optical Nn and Nm lines but 
overestimates significantly the Nii/16482 line, which leads to a 
larger uncertainty (0.3dex) on the resulting abundance (a mass 
fraction - 7 x 10" 3 ). Note that the [Nn] nebular lines were 
not employed in this determination. Upper limits to carbon and 
oxygen abundances of 0.07xsolar and 0.20xsolar were found, 
indicating a depletion in both elements. 

Our results agree very well with those found by Bresolin et 
al. ( 120021 for the WNLh star B 1 6 in NGC 300 and encouragingly 
- and as might be anticipated - land between the stellar parame- 
ters derived by Maryeva & Abolmasov (1201 11 1 for the minimum 
and maximum brightness phases (thus giving us confidence in 
our treatment of potential nebular contamination). Our derived 
E(B - V) — 0.06 is also consistent with the results of their anal- 
ysis and yields a deredenned V = 17.62. Interestingly, our de- 
rived Si abundance, which is roughly a factor of three lower than 
their value, agrees very well with the abundance of the other 
a-elements obtained from nebular abundances at that galacto- 
centric distance (Bresolin 2011). The nitrogen enrichment and 
the upper limits on the C and O abundances are consistent with 
an WNL evolutionary phase, characterized by an LBV-like be- 
haviour. 

Multi-epoch, quantitative analyses of the LBVs AG Car 
(Groh et al.l2009l . HD 5980 (Georgiev et al. I20TTT) and S Dor 
(Lamers 1995) have been made during their characteristic photo- 
metric excursions, and it is interesting to compare the properties 
of Romano's star to these objects. As with AG Car and HD 5980, 
the mass loss rate is found to increase and wind velocity decrease 
with decreasing temperature (Lamers 1995 does not report these 
parameters for S Dor). However, the bolometric luminosities of 
AG Car, HD 5980 and S Dor are found to decrease with the 
simultaneous reduction in temperature and increase in stellar ra- 
dius, whereas the opposite appears to be the case for Romano's 
star (Mareya & Abolmasov l201 11 1. 

It is thought that the reduction in bolometric luminosity for 
the first three objects results from the additional energy required 
to support the expansion of the outer layers of the stars; if cor- 
rect this would imply that the overal energy budget of Romano's 
star increased during its excursion. In this regard we highlight 
that the behaviour of the Galactic LBV AFGL2298 also implies 
an increase in the rate of energy production during its 'erup- 
tion', although in this case the growth in stellar radius occured 
at -constant temperature (Clark et al. 2009). If the results of 
these studies are replicated for these and other stars, the clean de- 
marcation between luminosity conserving excursions and non- 
luminosity conserving eruptions in LBVs will break down, de- 
spite the superficial similarities in their spectral and photometric 
evolution (cf. AG Car, S Dor and Romano's star); consequently 
raising the question as to whether a single physical mechanism 
underlies their divergent behaviour. 

5.3.2. M33 Var C 

M33 Var C has been the subject of considerable photometric and 
spectroscopic attention since its discovery (Hubble & Sandage 
119531) . albeit of irregular temporal coverage. Burggraf et al. 
(120111) present the most complete long term (-1900-2010) B- 
band optical light curve available, with observations reported in 
Humphreys et al. (Q988 ) and Szeifert et al. ( 1996) extending the 
wavelength coverage to the V- and R -bands. These data reveal 
two prominent photometric maxima (B ~ 14 - 15, V ~ 15.2) 
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between 1945-50 and 1984-8. Outside of these events, the star 
was typically observed with B ~ 16 - 17.5 (V ~ 16 - 17) al- 
though for at least 2 decades prior to its initial brightening in 
1940 it appeared to be in an extended faint phase (B ~ 17.5-18, 
V > 17). As can be inferred from the magnitude ranges quoted, 
significant short term variability also appears to be present (e.g. 
Fig. 4) in addition to the uncertainties associated with the histor- 
ical observations. Nevertheless M33 Var C currently appears to 
be in an extended bright phase which commenced in -1940 and 
upon which are superimposed two further outbursts separated by 
-40 yr (Burggraf et al. 1201 j} . The somewhat sparse post-1990 
photometric data presented by these authors appear to indicate a 
plateau state (B ~ 16.4) from -1990-2005 after which M33 Var 
C faded by ~0.8mag. 

Spectroscopic observations are likewise uneven. During both 
visual maxima M33 Var C demonstrated an F-supergiant plus 
Balmer emission line spectrum analagous to S Dor in its bright 
phase (Hubble & Sandage [19531 Humphreys et al. 1988). 
Outside of these events, the limited spectroscopy appears to in- 
dicate an earlier B spectral type dominated by Fe 11 and Balmer 
line emission - cf. the most recent high S/N and resolution spec- 
trum of Szeifert et al. (Q996), dating from 1992. More recently, 
Viotti et al. (2006) suggest a similar appearance in 2002 while, 
somewhat surprisingly, Burggraf et al. (1201 II) claim - but do not 
show - a much later spectral type coinciding with the current 
photometric fade. 

The spectroscopic and photometric data presented here sug- 
gest a hitherto unreported excursion analagous to the 1940s and 
1980s events. Commencing in -2001 July the optical peak was 
reached a year later and was maintained for a further year, dur- 
ing which M33 Var C demonstrated the F-supergiant emission 
line spectrum that characterised previous outbursts. If our in- 
terpretation of the mid-IR emission is correct (Sect. 4.2) this 
maximum light plateau endured until 2005 February. The lim- 
ited sampling of the B-band lightcurve of Burggraf et al. d201 11 1 
meant that they completely missed this event; unfortunately it 
casts doubt on their conclusion that the bright phases of M33 
Var C occur on a -40 yr period. Moreover, these data appear 
to demonstrate that spectral variability in M33 Var C may occur 
over relatively short timescales, with an F-hypergiant spectrum 
observed in 2003 (this work), B emision line in 2004 (Viotti et 
al. l2"0"06l late F in -2008 (Burggraf et aL l^OTlT) and B emission 
in 2010. For comparison recent complete outburst cycles (hot to 
cool phase and back again) of AG Car and HR Car occured over 
-10 and ~8 years respectively (Stahl et al. 120011 Szeifert et al. 
120031 . 



70, which require a dense and ionised medium to be driven into 
emission. Such a morphology is rare, being demonstrated by 
both M33 Var B and C (Szeifert et al. [l9"96l and this work) during 
their cool photometric maxima and the Galactic yellow hyper- 
giant (YHG) IRC+ 10 420 (Fig. 9), a star thought to be transiting 
rapidly to hotter temperatures as it executes a red loop across the 
HR diagram. 

However, the similarities to IRC+10 420 appear to end here. 
Prior to the 1990s IRC +10 420 underwent a long term photo- 
metric brightening (Patel et al. [2008 ) and while B324 also ap- 
pears photometrically variable (eg. Ma07) the limited data are 
consistent with a Cygni variability. Moreover IRC +10 420 is 
associated with significant dusty circumstellar ejecta of recent 
origin leading to a substantial IR excess (Blocker et al. 1999); 
in contrast, despite its intrinsic luminosity the mid-IR colours of 
B324 do not appear consistent with other massive stars experi- 
encing, or having recently undergone, heavy mass loss (Fig. 6Q 
Finally, B324 appears significantly more luminous than IRC +10 
420. Adopting V = 14.86 and (B - V) = 0.43 (Massey et al. 
120061) and the intrinsic colours of a F5 la star ((B - V) = -0.15; 
Fitzgerald [19701 leads to E(B - V) = 0.28 and A v = 0.87 (as- 
suming Ay - 3.1E(B - V) and subject to previous caveats re- 
garding the uncertainty in dereddening individual stars). Then 
assuming a negligible bolometric correction and /v=24.92 we 
find log(L»/L ) ~ 6.3, placing it significantly in excess of the 
empirical Humphreys-Davidson limit. 

Such a luminosity appears to be incompatable with a star 
evolving through a YHG phase via a red loop across the HR 
diagram, suggesting instead that an identification with a cool 
phase LBV outburst may be more appropriate (cf. M33 Var B 
and C). A consequence of this is that since there is currently no 
observational evidence for long term secular variability, B324 
appears to have remained in its current extended 'outburst' phase 
for at least twenty years. However, decades long outbursts have 
been observed in both the cool hypergiant M33 Var A (>35 yr; 
Humphreys et al. l2006l and R127 (-30 yr ; Walborn et al. 2008 ), 
while M33 Var C has yet to return to its pre-1940 quiescent 
state, suggesting that such behaviour is not without precedent. 
Intriguingly, the BHG [HS80] 110A is also located above the 
H P limit (log(L/L )~6.5±O.l and T eS - 21+lkK; Urbaneja et 
al. 120051 ) and shows no evidence for LBV-like variability over 
the past > 17yr (Sect. 3.3). Likewise, the Galactic BHG Cyg 
OB2 #12, which is also found to violate the HD limit, has re- 
mained stable for over a century (Clark et al. 120121 ). 



5.3.3. [HS80] B324 

Visually one of the brightest stars within M33, B324 is discussed 
at length by Monteverde et al. d 19961 ) who highlight both its late, 
possibly variable spectral type (F0-F5 Ia + in 1994 versus A5 Iae 
prior to 1990; Humphreys et al. 119901 ) and LPV in the Balmer 
series. P Cygni lines are observed in the Fell multiplets 40 and 



17 We suspect the IR emission of B324 arises from a combination 
of photosphere+wind, although interpreting the properties of such cool 
stars is not possible utilising current model atmosphere codes such as 
CMFGEN. For example, while we were able to reproduce the photo- 
spheric absorption spectrum of the LBV Wd 1-243 in the cool phase we 
were unable to simultaneously reproduce the emission lines originating 
in the stellar wind (Ritchie et al. l2009> . 
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Fig. 9. Comparison of the 1996 spectrum of the emission line F 
hypergiant B324 (solid line; Monteverde et al. 1996 ) against the 
Galactic YHG IRC +10 420 degraded to a comparable spectral 
resolution (red dotted line; Oudmaijer et al. 1998). Similar Fen 
P Cygni emission lines are also observed in the 2003 spectrum 
ofM33 Vai-C. 



6. Global considerations and concluding remarks 

Finally, what may we infer from the current population of can- 
didate LBVs within M33? Together with the photometric selec- 
tion criteria employed by Ma07, multiepoch observations em- 
ploying wide field imagers and multiobject spectrographs on 4- 
m class telescopes appear ideally suited to their study, as illus- 
trated by the identification of a third LBV excursion of M33 Var 
C between 2001-5. Combining existing spectroscopic and pho- 
tometric datasets reveals a number of new LBV candidates in 
addition to the well studied historic examples (e.g. Szeifert et 
al. 119961 Viotti et al. 120061 Ma07). The dramatic spectroscopic 
variability of the P Cygni star J013339.52+304540.5 and the 
WNL star J013332.64+304127.2 (Ma07) makes a compelling 
case for both objects. Pronounced long term photometric vari- 
ability is present in the BHG J013429. 64+303732.1 and the iron 
star J013459.47+303701.0. 

Moreover, low level optical photometric and spectroscopic 
variability appears ubiquitous amongst our sample, with LPV 
in the wind lines present in examples of all the (modified) 
classifications introduced by Ma07 and in the photospheric 
lines of the BHGs (particularly the Mg n/He i diagnostic lines). 
Unfortunately, our current temporal sampling is insufficent to 
distinguish between an origin of the former behaviour in short 
lived stochastic or long lived LBV driven changes in the 
properties of the wind, while the latter could arise in pho- 
tospheric pulsations, secular evolution of stellar temperature 
and/or changes in the structure of the wind/photospheric bound- 
ary. Nevertheless, optical photometry strongly suggests that 
many of the stars are a Cygni variables and hence likely pul- 
sationally unstable. 

We find that the spectroscopic data are also of sufficient qual- 
ity to permit non-LTE model atmosphere analysis and we present 
a proof-of-concept study of Romano's star. The successful anal- 
ysis of this star demonstrates the efficacy of our observational 
approach, but emphasises that high S/N and resolution data are 
mandatory to distinguish between Ofpe/WNL and P Cygni clas- 
sifications and to determine the temperature of the BHGs (e.g. 



via the Si n:Si m ratio). Moreover, more extensive spectral obser- 
vations would be required to ascertain the full elemental compo- 
sition of such stars as well as to drive down the uncertainties on 
parameters such as the wind clumping factor. 

As expected from the spectroscopic data - being intermediate 
between those obtained in photometric minima and maxima - we 
found that the physical properties of Romano's star in 2010 im- 
ply a smooth progression between these extremes. We note that 
the increase in bolometric luminosity with increasing radius and 
decreasing temperature is in the opposite sense to that demon- 
strated by AG Car, S Dor and HD 5980, in which bolometric 
luminosity is observed to decrease. In combination with analy- 
sis of AFGL2298 this suggests that the simple division of LBV 
behaviour into bolometric luminosity conserving excursions and 
non-conserving eruptions is overly simplistic. 

As suspected by Ma07, the subset of iron stars appear to be 
rather heterogeneous containing, in addition to strong LBV can- 
didates, a number of objects that can formally be classified as 
B[e] stars on the basis of their optical and mid-IR properties. 
Their nature is currently uncertain, but the lower luminosity sub- 
set may be related to the class of dusty, highly obscured objects 
identified by Thompson et al. (2009 ) as potential progenitors of 
transients such as SN 2008S, albeit suffering less self extinction. 
These authors find that approximately 25 of the 45 objects with 
similar mid-IR colours and magnitudes to the low luminosity 
B[e] stars (M [4 . 5] < -11.5 and 0.5<[3.6]-[4.5]<1.5) are like- 
wise optically faint or invisible, suggesting that they comprise 
an important component of the population of the most IR lumi- 
nous stellar objects within M33. 

Of the remaining iron stars, N93351, N125093, 
J013350.12+304126.6 and J013442.14+303216.0 possess 
optical spectra consistent with confirmed LBVs and mid-IR 
properties (colours and magnitudes) approaching known mas- 
sive evolved stars with extensive recent or ongoing mass loss; 
a connection with the recently discovered 'Object X' - the 
brightest stellar mid-IR source detected within M33 (Khan et al. 
1201 II ) - would be unsurprising. No stars appear to be as mid-IR 
luminous as rj Carinae currently is, despite over a century 
elapsing since its eruption. While this might be considered 
evidence of the rarity of comparable events, such a conclusion 
is hampered by observational bias - one might expect highly 
luminous stars to be located in active star forming regions that 
may themselves be IR bright leading to a lack of detection 
due to source confusion. Five further stars have 8^/m fluxes 
that are suggestive of the presence of dusty ejecta. While the 
bona fide LBVs M33 Var C, Romano's star and spectroscopic 
candidates J013339.52+304540.5 and J013332.64+304127.2 
lack near-mid IR excesses we may not at present exclude the 
presence of ejection nebulae associated with these stars, since 
emission from such dusty structures surrounding Galactic LBVs 
typically peaks at longer wavelengths. 

Adopting the modified classification scheme of Ma07 and 
excluding the non-LBV subset of BHGs and iron stars we in- 
fer a progression of Ofpe/WNL — > P Cygni — > iron star/BHG 
— > cool hypergiant for the transition from high to low temper- 
atures. Due to the empirical luminosity dependance of the am- 
plitude of LBV excursions, lower luminosity examples may not 
reach the Ofpe/WNL and P Cygni phases, while some well stud- 
ied LBVs such as AG Car have yet to be observed in a cool 
hypergiant phase (Stahl et al. 2001). Conversely, Romano's star 
and the SMC system HD 5980 (Koenigsberger et al. 120101 and 
refs. therein) suggest that LBV-like excursions may also extend 
to somewhat earlier WN subtypes than employed in this scheme 
(<WN8), while some cool hypergiants such as M33 Var A and 
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p Cas are instead observed to transit to cooler temperatures in 
outburst. 

In conclusion, such galaxy wide multiepoch observations 
employing large sample sizes appear well suited to the deter- 
mination of the duty cycle and physical nature of the outbursts 
and eruptions of LBVs and related objects, as well as the length 
of the LBV phase (via comparison of number counts to main 
sequence, Wolf-Rayet and cool super-/hypergiant populations). 
However, several caveats must be considered. Firstly, any such 
survey should ideally extend to earlier and later spectral types 
than included in the classification scheme of Ma07 in order to 
sample the full range of excursions undertaken by luminous 
evolved stars. Moreover, given the apparent LBV excursion of 
J013459.47+303701.0 (with V ~ 18.4) any study must also 
to extend to comparatively faint candiates which are currently 
poorly sampled. 

Combined with historical data, these observations support 
the finding that some excursions/outbursts of luminous evolved 
stars may last several decades. Both J013350. 12+304126.6 and 
B324 have spectra that resemble those of LBVs in outburst but 
neither have shown evidence of secular variability over the last 
~20 years. Likewise M33 Var C has yet to return to its pre-1940s 
quiescent state and consequently it could be argued to have been 
in 'outburst' for ~70 yr. Indeed, it appears possible that current 
assumptions as to the typical timescale of LBV excursions and 
eruptions have been biased to shorter durations by the length of 
the typical astrophysical career (as well as a small sample size). 
If such an assertion is correct, a considerable long term obser- 
vational investment will need to be made in order to understand 
such behaviour, and its relation to the more rapid (<yr) 'out- 
bursts' and pulsations that are exhibited by cooler transitional 
objects such as YHGs and RSGs (e.g. Clark et al. l20T0l Finally, 
despite their apparent stability, both B324 and the BHG [HS80] 
11 OA are found to exhibit combinations of temperature and lu- 
minosity which place them in excess of the HD-limit, joining 
a handful of other stars such as Cyg OB2 #12 that also appear 
to violate it; intriguingly the latter object also appears to have 
remained stable for over a century (Clark et al. I20121 >. 
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Appendix A: Presentation of spectroscopic data 
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Fig. A.l. Optical spectra of iron stars. A spectrum of the LMC LBV R127 obtained in during outburst in 1999 July is overplotted 
(red dashed line) on the spectrum of J013350. 12+304126.6 for comparison, after having the resolution degraded to match. 
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Fig. A.2. Optical spectra of BHGs. 
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Fig. A.3. Optical spectra of candidate P Cygni LB Vs. A high resolution spectrum of P Cygni itself is overplotted (red dashed line) 
on the spectrum of J013416. 07+303642.1 after having the resolution degraded to match. 
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Fig. A.4. Optical spectra of candidate Ofpe/WNVLs. A spectrum of P Cygni is overplotted (red dashed line) on the spectrum of 
Romano's star for comparison, after having the resolution degraded to match. 



J. S. Clark et al.: On the nature of candidate luminous blue variables in M33 



23 



Mi 



Var B 



Nil 



"ID QJ D CD 

XX XXX 



a 

M 

33 ~ 

I D 

33 



J013416. 10 + 303344. 9 



Fell - [FellJ 




i 



Var C 



I 




4 000 



4200 



4400 4600 

Wavelength (Ang) 



4800 



5000 



5200 



o 



J013416. 10 + 303344.9 



Var C 



r; 



173 



0) 

1 



6000 



6500 

Wavelength (Ang) 



7000 



Fig. A.5. Optical spectra of (candidate) LB Vs. Note that the spectrum of M33 Var B dates from 2003, and the spectrum of P Cygni 
has been overplotted (red dashed line) for comparison, after having the resolution degraded to match. 



